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Quantum systems subject to time periodic fields of finite amplitude A have conventionally been 
handled either by low order perturbation theory, for A not too large, or by exact diagonalization 
within a finite basis of N states. An adiabatic limit, as A is switched on arbitrarily slowly, has been 
assumed. But the validity of these procedures seems questionable in view of the fact that, as N ^ oo, 
' the quasienergy spectrum becomes dense, and numerical calculations show an increasing number 

of weakly avoided crossings (related in perturbation theory to high order resonances). This paper 
deals with the highly non-trivial behavior of the solutions in this limit. The Floquet states, and the 
associated quasienergies, become highly irregular functions of the amplitude A. The mathematical 
radii of convergence of perturbation theory in A approach zero. There is no adiabatic limit of the wave 
functions when A is turned on arbitrarily slowly. However, the quasienergy becomes independent of 
time in this limit. We introduce a modification of the adiabatic theorem. We explain why, in spite 
of the pervasive pathologies of the Floquet states in the limit oo, the conventional approaches 

are appropriate in almost all physically interesting situations. 



> 

(N 
OO 



s 

I 



42. 50. Hz, 42.65.Vh, 03.65.-w, 05.45. +b 



O ■ I- INTRODUCTION 



Physical systems subject to finite time periodic perturbations of amplitude A and period T = It: jio have been 
extensively studied by making use of the Floquet theorem. This theorem, a consequence of the discrete time 
translation symmetry of the Hamiltonian, states that there is a complete set of quasiperiodic solutions of the time 
dependent Schrodinger equation which, when i — > i -I- T, are simply multiplied by a phase factor exp(— ze„r), where 
'■^ e„ is called the "quasienergy". This phase factor defines e„ only modulo lo, and so the quasienergy may always be 
^ \ taken to lie in the strip < e„ < lo. 

O ■ For such systems the quasienergies e„ are of comparable interest to the energy levels of time independent systems. 

\ One would expect to use the subscripts n as unambiguous labels of the time dependent states of a periodically driven 
^ . system as the magnitude A of the perturbation is switched on adiabatically, as in the case of time independent quantum 
systems. Indeed, there is a substantial literature studying Floquet systems along these lines, usually employing 
numerical methods, in a finite set of basis states ||^. Of special interest are the "avoided crossings" (ac's) — regions 
in the A — e plane where two quasienergies approach each other as a function of A and (except for special symmetries) 
avoid crossing one another. 

But there is a difference in principle when a complete infinite set of basis states is included. For a spatially confined 
system, with an infinite number of discrete energy levels for A = 0, there is also an infinite number of quasienergies 
for A > 0, and the spectrum fills the fundamental strip densely. In fact, typically it is a dense point spectrum 
(For example, in the special case of a particle in a one dimensional square well and vanishing A, Weyl |^ has 
shown that the energy eigenvalues fill the fundamental strip densely and uniformly) . As the number of basis states 
becomes infinite there is a weakly avoided crossing (hereafter "wac") near each point in the A — e strip. This leads 
to qualitatively new issues: Do the Floquet state solutions for a given A converge to well defined limits? Are the 
Floquet states and quasienergies well behaved functions of A? And is there a well defined limiting path which the 
system follows as A is switched on arbitrarily slowly — i.e., does an adiabatic limit exist? We have examined these 
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questions in a variety of approximate ways and have arrived at a coherent picture, though generally wo do not have 
conclusive mathematical proofs. We find that the radii of convergence of power series expansions in A, starting 
from the unperturbed eigenstates, are zero. Floquet states and their quasiencrgies are discontinuous functions of A 
everywhere; there is no adiabatic limit in the usual sense (though we will propose a useful weakencxl inodification of 
the adiabatic theorem). At the same time, in the limit of slow switching on of A the quasienergy remains arbitrarily 
close to its initial (A = 0) value. We will explain the consistency of these results with the well established success of 
standard time dependent perturbation theory and of the adiabatic theorem. 

We note that these features will not be seen directly in any numerical study as N becomes large. The effects become 
increasingly weak very rapidly as the basis size is increased. Although, for a given interval in A, inclusion of very high 
lying levels does have a major impact on Floquet states and quasienergies, it is only over an increasingly smaller range 
of A, and this becomes at some stage unobservable on the scale of numerical accuracy available to the computer. 

In Section II we review elementary Floquet theory in a finite basis. The problems arising for an infinite basis 
are discussed in Section III. In Section IV we analyze the convergence of the Floquet states as the size of the basis 
becomes infinite. Section V is devoted to the dependence of the states and quasienergies on A, including questions of 
labelling Floquet states, and of the existence of an adiabatic limit, when N oo. The relationship to finite order 
time dependent perturbation theory is discussed in Section VI. Our conclusions are summarized in Section VII. 



II. FLOQUET THEORY IN A FINITE BASIS 

In this section we review the elements of time periodic Hamiltonian systems, 

H{t) = Hq + XHi{t) with Hi{t + T) = Hi{t), (2.1) 

when approximated by a finite matrix, H^{t), in the representation of the A'' lowest states of the time independent 

Hamiltonian Hq. In the next section we will discuss the considerable problems that occur in the limit N oo. 

Because of the discrete time translational symmetry of H^{t), there is a complete set of solutions f^it) of the time 
dependent Schrodinger equation (with h= 1), 

-•^^"^ -il^(0<p^, (2.2) 



which are of the Floquet form. 



with a time periodic part 



dt 



^^{t)=eM-<t)u^it), (2.3) 



<(t + r)=<(i); (n=l,2,...,Ar), (2.4) 



and "quasienergies" which may be taken to lie in the interval [0, w), with u = 2-k/T. For given A one may label 
the Floquet states in order of increasing quasienergy, the eigenvalue which characterizes discrete time translational 
symmetry. 

Figure 1 shows such a quasienergy spectrum (for finite N) as a function of A for the example of a free particle in a 
one-dimensional box [x G [—a, a]) with harmonic driving, 

iJ(t) =//2rn + Asin(^)cos(wt). (2.5) 

The quasienergies are continuous functions of A that do not cross, but show avoided crossings (ac's) provided there 
are no symmetries which allow actual crossings. These ac's are abundant in the spectra of time periodic systems and 
are of central interest in this work. As the basis size is increased, the newly introduced avoided crossings tend to 
become rapidly weaker (smaller gaps at the crossing and smaller range of A over which their effects are siibstantial). 
They therefore tend to become isolated from one another, and it is useful and meaningful to consider their effects 
individually, as we do in the following. 

As two quasienergy lines pass an isolated weakly avoided crossing (wac) at A = Ai , the corresponding Floquet states 
rapidly interchange their forms. At the point A = Ai they are very nearly linear combinations of the two Floquet 
functions just outside the region of the ac, with amplitudes of equal magnitude. As |A — Ai| grows, this mixing 
decreases rapidly, as follows. The functions are mixed by no more than a given small relative weight 5, provided that 
|A- All > w{5)/2, where 
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^(-^^ = Ai r (2.6) 

0|cri - (72 I 

Here Ae is the quasienergy splitting at the ac, and cri and (T2 are the slopes de/dX of the quasienergies at the crossing 
point in the absence of the terms connecting these states (see Fig. 1). We neglect the weak influence of all other states 
near this wac. 

Using standard Floquct state perturbation theory |^ in A, one finds that the second order term describes very well 
the overall behaviour of the quasienergies for small A. However, near A = Ai nearly degenerate perturbation theory 
is needed, giving 
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ei,2(A) « - [((71 + (72)(A - Ai) ± V('^i-^2)2(A- Ai)2 + (Ae)2j . (2.7) 

It is helpful to consider complex values of A. Then the time evolution is no longer unitary, and quasienergies are 
complex. They may be considered as the N values of a single TV-valued analytic function, e(A), with N Riemann 
sheets connected at complex branch points 0]. As is the case for real A, the function e(A) is defined only modulo w, 



and we choose always < Re e(A) < uj. From (2.7) we see that wac's for real A are manifested as branch points near 
the real axis of e(A), at 

A = Ai ± , A . (2.8) 

(see Fig. 2). When A passes along the real A axis through Ai, as discussed above, there is a rapid change in the spatial 
part of the two Floquet states, which are approximately interchange d. In contrast, on a path starting from the real 



axis, going in a loop in the complex plane around the branch point (2.8) and back to the real axis, each of the two 
Floquet states returns approximately to its original spatial dependence. Perturbation expansions in A have finite radii 
of convergence. At A = the eigenstate cpj with energy Ek defines the quasienergy ej = Ek (mod lo) on the j-th 
Riemann sheet of e(A) (the indices k and j are unequal, in general, because j labels increasing values of e in [0,a;), 
whereas k labels increasing energy values E over [0, 00) ). Then the radius of convergence Acj- for ipj and tj is the 
magnitude of A at the branch point nearest the origin on that j-th sheet. 

In the laboratory the perturbation, A-ffi, is commonly turned on slowly. This can be characterized, as usual, by a 
switching factor e** in the interaction (s > 0): 

H{t) ^Ho + Xe^'Hiit) where Hi{t + T)= Hi{t), (2.9) 

with the initial condition 

^fi-^) = ify (2.10) 
We are interested in the wave function V'J'^(^) at a specified time, say t = 0. This system is no longer periodic in 



time, and the solutions of the time dependent Schrodinger equation corresponding to (2.9) and ( 2.10| ) are therefore no 
longer Floquet functions. However, if the turn-on rate s is slow compared to the driving frequency uj, it is useful to 
describe the solutions at time t in the basis of the Floquet functions (2.3) at the corresponding value of \{t) — Ae**. 
Two results are known ^ : 

i) Every ac can be characterized by a rate 

|(7l - (72|Al 

For large turn-on parameters s 3> C (provided that s <^ lo) the solution will follow closely the initial Floquet state, 
as if there were no ac. This is the so-called Q Landau-Zener transition through the ac. On the other hand, in the 
adiabatic limit for this ac, s ^ ^, it will closely follow the Floquet state which passes the ac continuously. For s w ^ 
the solution, after A(t) has passed Ai, will be a superposition of the two Floquet states involved in the ac. 

ii) The adiabatic theorem: Apart from an overall phase factor, the final state ^'^''''(0) converges in the limit s 0, 
and the limit is the Floquet state of the periodic system that corresponds to the quasienergy found by following the 
quasienergy ej of the initial state ipj as a continuous function of A. 

These properties are well established for Floquet systems with a finite basis. In the next section we show that the 
limit N ^ 00 \s highly pathological. 
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III. THE LIMIT N ^oo 

The general difficulties raised by an infinite basis are outlined here; a detailed discussion is given in the following 
two sections. 

For definiteness we will restrict the discussion to finite one-dimensional systems with x G [—a, a] and 

H{t) = p^/2m + U{x) + XV{x) cos(cji), (3.1) 

where U{x) and V{x) are analytic and bounded. Re A will be kept uniformly bounded, Re A < Ai, as iV becomes 
infinite. For A = the conventional energy spectrum Ek is discrete, and for sufficiently high eigenvalues exhibits 
increasing spacings between successive levels (approximately proportional to the level index, as for the case 
U{x) = 0). We will make repeated essential use of this feature. 

In the limit N —^ oo one expects that for every A the N quasienergies within the finite interval [0,0;) will form a 
dense spectrum. In fact, Weyl showed [Q that for t/(x) = A = 0, so that Ek — (l/2m)(7r/2a)^fc^, the spectrum is 
uniformly dense, provided only that uj/Ei is irrational. This is the generic situation Q for U{x), A 7^ 0. Moreover, 
Howland ||] has shown that, for Hamiltonians of the form (^^), and for most values of the coupling strength A, the 
quasienergies have a dense point spectrum. 

What happens to the ac's in the limit N —^ oo7 We find a very simple picture in the complex A plane: ac's correspond 
to branch points, as in Fig. 2. With increasing TV we add higher lying states of the unperturbed Hamiltonian to the 
basis and find that the gaps of newly introduced ac's, as well as the imaginary part of the corresponding branch points 
(see the Appendix), decrease faster than any power law with N . Therefore the real A axis is a line of accumulation for 
the branch points. Moreover, the branch points with Im A > A2 are finite in number and tend to well-defined limits 
as iV — > 00, for any A2 > 0. The problems with the limit N ^ 00 are restricted to the immediate neighborhood of 
the real A axis. 

For real A one expects the ac's to be dense in the A — e plane. In fact, along each quasienergy curve e(A) obtained 
within a finite basis approximation one expects a dense set of ac's when N 00. This has a number of consequences: 

i) As there will be ac's on each quasienergy curve for arbitrarily small values of A, the radii of convergence Xcj of 
the perturbation theory in A for quasienergies and Floquet states all shrink to zero. 

ii) The increase of the basis size from an initial value N will introduce ac's arbitrarily close to any point on a 
given quasienergy curve defined with the initial basis of size N. Each of these ac's will make large changes in the 
corresponding Floquet states over a finite A-interval, strongly admixing and interchanging pairs of states, as described 
above, and it is therefore by no means clear whether or not Floquet states converge as A'^ — > 00. 

iii) For any finite N one can label the Floquet states iy9^'^ and quasienergies e^'^ in such a way that they are 
continuous functions of A (if, as usual, we identify e with e + lo). If we increase N, we have to rearrange labels for 
every new ac, and it is not clear if there exists a meaningful labelling which tends to a well defined limit as N —^ 00. 

iv) Let the periodic field be switched on over the time interval —00 < t < 0, with X{t) = e**A. For any finite N 
there is a well defined adiabatic limit, as s — > 0, of the state (up to an overall phase factor) and of the quasienergy at 



the final time, t = 0. Let ^min be the smallest of the rate parameters ^, defined by (2.11), characterizing the relevant 
ac's in [0, A]. Then, for s <C ^min the solution of the Schrodinger Equation ( |2.9| ) will simply follow the Floquet state 
corresponding to a continuous quasienergy curve. In the limit A^ — > c», however, there will be ac's with arbitrarily 
small parameters ^. Therefore, for smaller s more and more of these weak ac's will be passed adiabatically, rather 
than undergo Landau-Zener transitions, leading to completely different final states tpl,'^{t). Thus an adiabatic limit 
as s ^ cannot be expected. 



IV. CONVERGENCE AND NON-CONVERGENCE OF FLOQUET STATES 



For systems of type ( B.l ) we shall present strong arguments that, for a set of full measure in A in an A-independent 
interval < A < A, the Floquet states converge in the limit A^ 00, even though the ac's are dense in the A — e 
plane. This does not imply, however, that the limit is a continuous function of A. For we will also argue that the ac's 
give rise to infinitely many A's, of measure zero but dense in any A-interval, for which none of the Floquet states of a 
finite basis converge. This can be stated more precisely as follows: 

Proposition I (Convergence): For any interval [Aa,Af,], any 6 > 0, and any 77 > 0, there is an integer 
M((5, 77, Aa, Af,) with the following properties: For any M' > M and each m = l,2,...,Af there exists a label 
m' G {1, 2, . . . , Af'} and an overall phase factor e*" such that 

|^i^'^x,i)-e-^5''(x,t)|<5 (4.1) 
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and 



|expH(e^'''-e™''')r]-l|<2<5, (4.2) 



for all X, for all t £ [0, T], and for all A G [Aq, Af,], except for a subset of A of measure < rj. 

Proposition II (Non-convergence): For any interval [Aq, Ab], any S > 0, any M, any label m £ {1,2,..., M}, 
sufficiently large M' > M and any label m' e {1, 2, . . . , M'}, 



1 

2a 



< 5, (4.3) 



for a set of A's dense in [Aa, Af,] and all t. That is, for this dense set of A's the eigenfunctions within the smaller basis 
have arbitrarily small overlap with any of those corresponding to the larger basis; they do not converge to a limit. 

Proposition I states that as N grows, there is an increasing measure of A where a finite basis calculation gives 
the Floquet states and quasienergies correctly within an arbitrarily small error. For many practical purposes this 
supports the use of a finite basis for describing a Floquet system. Nevertheless, it is important to realize, as stated 
in proposition II, that even for an arbitrarily large basis size N there are infinitely many A's, dense in any interval 
(albeit of total measure zero), where a given Floquet state (p^ (t) does not converge in the limit N — > oo. 

We can prove proposition I for any system with the follo wing model property which, we submit, captures for this 
purpose the essence of a real Floquet system of type ( |3.l| ): For an ac at A„ £ [Aq, Xb] we again define the interval 
A„ ± t«„((5„)/2, outside of which the admixture of the two unperturbed Floquet states changes by no more than (5„. 
Then we assume that we can choose a set {Sn}, with < ^ for any chosen 5 > 0, such that J2'^=i "^niSn) 

converges. This assumption appears to be satisfied in systems of type ( |3.lD , although we have no mathematical proof. 
Numerical calculations and analytical considerations (see Appendix) suggest strongly that the gaps Ae„ decrease 
faster than any negative power of n, due to the increasing spacings of the unperturbed energies Ej. Let us choose 
6n to dec rease relatively slowly, as a small power of n, say (5„ = (5/(2n^). Then the sum over Wn will converge (see 
Eq. (|2.6| )), as long as the difference |cri — (T2I in quasienergy slopes doesn't decrease too rapidly with n. Numerical 
experience suggests this to be the case. 

We make an argument based on the Borel-Cantelli lemma |llT| . Let us start, as usual, with an approximation to 
the system given by restriction to a finite number M of spatial basis functions. The new ac's introduced as this is 
increased to a complete, infinite, basis set are labelled from hm to 00. For a large enough choice of the initial basis size 
M, the partial sum J2'i^=nM ''^ni^n), which gives the measure of A's where a Floquet state might be affected by more 
than S when increasing the basis size from M to infinity, can be made smaller than any given ry. For all other values of 
A, those that are within none of the intervals A„ ± ■u;„((5„)/2, and which therefore constitute a s et o f measure at least 



A;, — Aa — r], the Floquet states are changed by no more than J2n=i ^- This explains Eq. (4.1) of proposition I. 
Equation (4.2) follows at once from Eq. (4.1). The quasienergies can be determined from e""^"^ = (p^\x,t + 



T) / ip^^ {x , t) . For Floquet states normalized by l/(2a) J^^dx\ip^^ {x,t)\'^ = 1, we can choose an x and t where 

|iy9^^(a;, i)| > 1. Th en can be determined from 93^^ {x,t), which differs from (p^^{x,t) by less than S (Eq. (4.1)), 
leading to Eq. (|0). 



Proposition II follows from the plausible, but unproved, assumption that any quasienergy line of a finite basis 
(Af ) calculation will show an ac within any given A-interval, if the basis size is increased sufficiently. (Failure of this 
assumption would imply that there do exist finite A-intervals in which a quasienergy line is never crossed as A/ — + cx3.) 
Using that assumption, we can argue straightforwardly: Within any given interval / = [A', A"] C [Aa, A;,], for a suitably 
large basis size M one will eventually find an ac. This crossing changes the ni-th Floquet state (p^^^ by more than 
some chosen amount, say 40% admixture of orthogonal basis states, over some finite A-interval. Within that interval 
the two Floquet states of the ac then have an overlap of less than 0.6 with (p^^ . We now select one of these states. 
With a further increase in the basis size, it will ultimately encounter an ac within the chosen interval and will be 
changed by more than 40% over some smaller but still finite A-interval. Within that interval for an even larger basis 
size also the second state will be changed by more than 40% due to an ac. This leads to four Floquet states within 
a finite A-interval, each with an overlap of less than (0.6)^ with (p^^^ . Repeating this argument sufficiently often, one 
finds a basis size Mi > M and an interval J„j C / where no Floquet state ip^^'^ has an overlap of more than a given 5 
with (y9^^ of the initial basis size M, leading to proposition II. As stated in proposition I, these A's, even though they 
are dense, have a measure which tends to as M — > 00. 



We can derive a stronger version of proposition II. Repeating the argument above for the n2-tli Floquet state 95, 



M 

"2 



of the initial basis size M, one finds a basis size M2 > Mi > M and an interval /„2 C /„i C / where no Floquet state 
93*^^ has an overlap of more than a given 6 with 93*^ nor with ip^^^ . This argument can be repeated for all Floquet 
states iy9*^ of the initial basis. Thus there is a dense set of A's in any A-interval, where none of the Floquet states of 
a finite basis converges [fH . 
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V. FUNCTIONAL DEPENDENCE ON PERTURBATION STRENGTH; ADIABATIC LIMIT 



A. Labelling of Floquet States 

Proposition II of the last section has an immediate consequence: A continuous labelling of the Floquet states y^'^ 
and quasienergies e^'^ as a function of A, possible for finite N, is no longer possible in the limit N ^ oo. We therefore 
propose a new way of labelling, useful at least for small A: We assign the label n to a Floquet state, if its overlap 
with the n-th eigenstate (pn{x, t) of Hq is larger than 50%, in the sense 

'^^ ^ / dx \^^{x,tr^^,{x,t)\ > 0.5. (5.1) 

This procedure will not always find a label for a Floquet state. This is obvious for large A, where none of the Floquet 
states resembles a low lying unperturbed state. But, even for small A, in the center region of an ac the overlap with 
an unperturbed state will be less than 50%. 

Since for small A the quasienergies are very flat as a function of A, and ac's will generically occur only between 
Floquet states related to states of Hq that are far apart in energy, their ac's are expected to have a width that 
decreases faster than any power law as A goes to zero (Appendix). We therefore expect that, on the interval [0, A[, 
the labelling works for a Cantor set of A values with finite measure less than A, and that this measure approaches A 
as A ^ 0. 



B. Adiabatic Turn On 

Even though perturbation theory for Floquet states and quasienergies in A has zero radius of convergence in the 
limit N — > oo, we find a simple, strict result for the perturbation expansion of a solution ?/;^'^(i) of the Schrodinger 
equation ( |2.9| ) , where the periodic driving is turned on from t = — cx) by the factor Ae^* (we remark again that V'fc (0 
is not a Floquet state). 

Theorem: The perturbation expansion of ip^''^{t) in A, for t < oo, has an infinite radius of convergence for any 
s > 0, i.e. is an entire function of the complex variable A. 

This can be proven by majorizing the perturbation expansion. With the wave function x{t) expressed in the 
interaction picture, 

x{t)=e^"ot^l'\t), (5.2) 

the Schrodinger equation (U) becomes idx{t) /dt = XW{t)x{t), with W{t) = e'^^^e^^Vix) cos{ujt)e-'"°* . The famiUar 
formally iterated solution is 



l + i-iX) I dtiW{ti) + {-iX)^ I dt2 r dtiW{t2)W{h) 

CO J —oo J —oo 

t rt3 rt2 



oo ^ — oo ^ — oo 



{~iXf / di3 / dh / dtiWm)W{t2)W{t,) 



X{t = -CX3) (5.3) 



Since e is unitary, it can easily be shown that for any normalized states / and g 

<f\W{tr.)W{t„^,)...W{h)W{h)\g> < Caxe^^*"+-+*^\ (5.4) 

where Knax is the maximum of |l^(a;)|. From this the n-th order term of the perturbation expansion of < f\xit) > 
can be majorized by (l/n!)(AT4iax/s)"e'** and thus the expansion ( ^.3[ ) converges for any s > 0. 

That is, the state which evolves with a given switching-on rate s is uniquely and well defined (in a finite or infinite 
basis), for an arbitrarily large final interaction strength, in spite of the convergence problems with Floquet states. 
There is, however, an anomaly: there is no well-defined adiabatic limit when s — > (see also Sec. III). 

Although an adiabatic limit in the usual sense does not exist, we now show that for sufficiently small A there is a 
(logarithmically) large window of turn-on parameters s < s < s, where the final state V'fe^(^) is almost independent 
of s. 

Proposition III: For any small S > and large 77 > there exists a A > and an interval [s, s] with 
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s/s > 77, 



(5.5) 



such that for all s' , s" G [s, s], all A < A, all x and all states fc, 

^-<5, (5.6) 



A 

where e'" is a (physically uninteresting) overall phase factor. 

That is, for any desired level of convergence (as defined by 5) we can find a range of s of arbitraily large relative 
size (arbitrarily large s/s = 77) by restricting A to a sufficiently small value. 

Consider any finite N and the Floquet state arising out of the unperturbed state (pk- The perturbation series ( |5.3| ) 
involves first order terms of the form 

x"*'^-^El? (5.7) 

^ Ek' - Ek±uj - IS 

k' 

From this it is clear that to achieve independence of s with an accuracy (5, s must be smaller than a value s given by 

Vk'.k 



-sxy: 

k' 



{Ek' ~Ek±oj) 



rfk' 



< 6. (5.8) 



For small enough A all higher order terms up to any finite order n may be neglected (assuming, as we do, that there 
is no exact resonance Ek' — Ek = imuj, m < n). 

However, for any given Aq, no matter how small, there will exist a sufficiently large basis size such that there are 
some (weak) ac's associated with any quasienergy in the A-interval [0, Aq]. Consider for any initial basis size A^o an 
arbitrary quasienergy curve efc(A), continuous over the interval < A < Aq. As the basis size is then increased to a 
sufficiently large size TV' > A'o, avoided crossings of efe(A) will be introduced within this interval. These ac's will be 
characterized by rates (see Eq. ( |2.1l| )) which we label ^1,^2, ■ • • , in order of increasing energy of the states at A = 
from which the crossing curves arise, which then assures that they are ordered with decreasing rates: ^i+i < As 
N ^ 00 and A(i) (= Ae**) grows from to A all these (infinitely many) weak ac's will be encountered, and each of 
them gives rise to an admixture of a new state into ip'^'^{x,t = 0) with amplitude (see [^) smaller than 

Sn = V^Cn/s- (5.9) 

To make the total variati on in tp^^'^{x,t — 0), as defined by X^n^"' smaller than a given S over the whole range of 
turn-on rates set by Eq. (^.5|), we use the following properties of the ^'s in the limit A — > 0: 

lim = (5.10) 

A — >0 



limCi+i/ei =0 

A — >0 



(5.11) 



Equation ( ^.10 ) follows from the fact, that for small enough A the largest a c of in [0, A] will be with an arbitrarily 
high-lying state k' of Hq and that the rate £ 1 of the ac, according to Eq. (2.11) (see also the Appendix), decreases 
faster with k' than any power law. Equation ( ^.11 ) is due to the fact, that the ratio ^j+i/^; is, in the most unfavorable 
case, due to ac's of neighbouring levels k' + 1 and k' of Hq with e^. The ratio of the corresponding quasienergy 
splittings decreases exponentially with (fc' -I- 1)^ — k'^ — 2k' + 1 and goes to zero in the limit A ^ 0. 

From these properties we conclude that for s larger than any s < s/rj the variation in tfj'J,'^{x, t — 0) due to infinitely 
many ac's can be made smaller than any given 6 for sufficiently small A. Therefore for small enough A one finds for 
all A < A almost adiabatic behaviour in a window [s, s] of turn-on parameters with s/s arbitrarily large. 



C. Conservation of Quasienergy 

Here we consider again the limit of turning on X(t) arbitrarily slowly from up to some arbitrary A. We have 
seen that the eigenfunction V'fc^(O) does not have a limit for s —> 0. However, we will argue here (but not prove 
mathematically) that, as s ^ 0, V'fe'^(O) is within arbitrary accuracy a linear combination of Floquet states of the 
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periodic Hamiltonian (Eq. ( |2.l|) ) which have quasi ener(7ies arbitrarily close to the initial value efc(A = 0) — i.e., the 
energy modulo uj. Thus in the limit s — > quasienergy is a conserved quantity: 

Proposition IV: For any 6 > and any A > there exists an sq > such that for all s < Sq 

r dx\^l'\x,t)^e'"'^i:l'^ix,t + T)\^ <S (5.12) 

J — a 

for all t < 0, and for all states k. This is true even though lims^Qipl.''^{x,t) does not exist for any finite time 
-oo <t<0. 

To make this plausible we will use a simplified geometrical picture, in which quasienergies are linear functions of 
A. Although all ac's are now represented by actual crossings, they will be traversed dynamically like real ac's with 
some finite rate parameters ^. The infinitely many quasienergy lines fall into two classes, according to whether the 
magnitudes of their slopes are smaller or larger than a specified critical value CTc = <5'/(2A). We will further assume 
in this model that the latter (large slope) class, for any S' > and any A > 0, has only a finite number K{S', A) of 
members. 

In the course of turning on the periodic driving from up to A with a small enough turn-on parameter s, the 
infinitely many quasienergy lines with slopes < CTc will, at most, cause a deviation in quasienergy by UcA {— S' /2) 
away from the initial value efe(A = 0). The additional changes in quasienergy from ac's with the finite number K{5' , A) 
of steeper quasienergy lines can also be restricted to be less than 6' /2 by a sufficiently small choice of turn-on rate 
s: Since the ac's are dense on any quasienergy line, for sufficiently small s one will be diverted from any of the 
steep quasienergy lines within any given small quasienergy range, which we choose to be S' /2K. The total change in 
quasienergy as A is increased from to its final value thus can be made smaller than d'/2 + K6'/{2K) — 5' for any 
given 5' . For 5' = VS/T and within this simplified model we have thus proved proposition IV. 

The same arguments should hold for the original Floquet problem, as including the quasienergy dependence on A 
at the ac's reduces their steepness and thus further reduces the spread in quasienergy from the above estimates. Also, 
the overall nonlinear dependence on A which follows, e.g., from second order perturbation theory, poses no problem 
for the argument. The assumption that the number K{5' , A) of slopes larger in magnitude than Uc is finite, remains 
reasonable, since an increase in the basis size introduces Floquet states originating from higher lying states of iJg, 
which show decreasing dependence on A. However, the problem of defining slopes at all (in the limit N oo), as 
there is no continuous labelling of the quasienergies as a function of A, will make a mathematical proof difficult. 

This proposition leads us to a better understanding of the nature of the state il^l'^iO) in the limit s ^ 0. As i^k^i^) 
is a linear combination of Floquet states with quasienergies closer and closer to £k , it changes constantly as 6 and s 
go to zero. Thus, while there cannot be an adiabatic limit for the wave function, quasienergy is conserved in the limit 
s ^ 0. 



VI. THE STATUS OF TRADITIONAL FINITE ORDER PERTURBATION THEORY 

Non-linear optics is a major field of science in which traditional, finite order perturbation theory in the applied 
electric field (usually to low order) successfully describes experiment. Here we shall show why this well-established 
theory is consistent with our considerations in spite of our conclusion that, strictly speaking the radii of convergence, 
Xcj , of perturbation theory vanish. 

We reiterate first that if the turn-on rate s — > 0, pcrturbative non- linear optics in fact fails. For small A this failure 
is due to near-resonances, Ek' — ~ into, generally with very high-lying excited states. (This is the reason why, for 
any finite basis size N , Ac is finite.) 

When the perturbation is turned on as in Eq. (^^), we have seen that, provided the turn-on rate s is small enough 
but exceeds a lower limit, s, then as A — > 0, the resultant state can be made arbitrarily close to the traditional first 
order pcrturbative solution. A similar result can be derived for the traditional pcrturbative solution up to any finite 
order. In typical laboratory situations we have seen that s is exponentially and unphysically small. For finite small A 
the perturbation expansion is asymptotically convergent. 

Finally we briefly mention the unavoidable effects of line broadening. The quantum system of interest is inevitably 
in contact with its environment, and there are interactions between the many particles that ordinarily constitute the 
quantum system of interest, so the individual particle states are lifetime broadened. We conjecture that if broadening 
is characterized by a finite width F, then a finite radius of convergence will be restored. 
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VII. CONCLUSION 



There have been two standard approaches to deaUng with the behavior of quantum systems subject to strong 
time periodic fields. One is the use of finite order perturbation theory (e.g., second or third order nonUnear optical 
susceptibilities), and the other the exact solution of the problem within a finite basis of states. But both of these 
approaches miss qualitative features of the exact mathematical solutions. 

We have shown by a set of "propositions" (as opposed to rigorous mathematical proofs) that for a large class of time 
periodic problems the structure of the exact states and the quasienergy spectrum is remarkably irregular. By "exact" 
we mean here that the complete infinite set of basis states is included. Interaction with the environment is neglected. 
We have considered the states and quasienergies as functions of the strength A of the time periodic potential, as the 
number N of basis states becomes infinite. We have found that in any interval Ao < A < A;,, although the states 
converge to a well defined limit as A'' ^ oo for a set of A with the full measure Ab — Xa of the interval, there is a set of 
A, of total measure zero, but dense within every finite interval, for which the states do not converge. As a result, in 
contrast to the situation for any finite N, it is impossible to label states and quasienergies continuously as a function 
of A. The familiar quasienergy "dispersion" curves as functions of A (as shown, e.g., in Fig. 1) become discontinuous 
everywhere. One consequence of these discontinuities is the absence of a true adiabatic limit; there is no unique final 
state to which the system tends as the periodic perturbation is switched on arbitrarily slowly. 

But these pathologies, including a radius of convergence of perturbation theory in A which approaches zero as 
iV ^ 00, do not show up under most physically realistic circumstances. In particular, we have explained the familiar 
and well established success of ordinary time dependent perturbation theory in terms of the modified adiabatic theorem 
and the typical smallness of the parameter s which enters that theorem in practice, as well as the successes of finite 
basis calculations. 
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APPENDIX: EXPONENTIAL DECREASE WITH BASIS SIZE N OF NEWLY INTRODUCED 

QUASIENERGY GAPS 



We label the eigenstates of the time independent Hamiltonian (A = 0) by an index j which increases with the 
unperturbed energy. We consider the solution of the full time dependent problem in the limited spatial basis of the 
first N such states. We demonstrate here that, for sufficiently small A, the new ac's introduced by the inclusion of the 
next basis state (labelled 1) are characterized by quasienergy gaps that are smaller than a bound which decreases 
exponentially with N. 

As in Sec. II, let us take as the "iV*'' level" Hamiltonian the representation of H{t) in the basis of the first N 
unperturbed states: 

N 

H^{t)^H,+ J2 \i){i\V{m{j\^Ho + V''{t), (Al) 
where the unperturbed Hamiltonian is 

oo 

Since has the same time periodicity as the full Hamiltonian, the corresponding time dependent Schrodinger 
equation has solutions of the standard Floquet form l\2.. 



^j:{x,t)=exp{-ie'^t)u'^{x,t), (A3) 

with k = 1,2, ... N, where the functions {x, t) are time periodic. Therefore, the solutions formed from the basis of 
the first N states are of the form 



N oo 

E E «f(j;")e— (A4) 

j—l n= — oo 



where normalization imposes the restriction 

N oo 
j — l n— — oo 

Now we include in ( |A1[ ) the next highest state + 1), whose unperturbed energy -Ejv+i can be written as 

K+i = e + Mnco, (A6) 

with e confined to the fundamental stripe, < e < lu. At some value of the coupling A the quasienergy of this state, 
e, may equal that of one of the solutions labelled k (1 < fc < N) in the basis of the first N states (see Eqs. (A3) and 
( A4)). The perturbation of the remaining potential, V — , turns that into an avoided crossing, with a gap given 
approximately by twice the corresponding matrix element, 

N 

A = aE(A^ + mj) K(j; + 1) + (j; Mat - 1)] . (A7) 

We now place strong limits on the size of the right hand side of this equation. Since the state label k and the basis 
size N will remain fixed, for simplicity of notation we will no longer write the subscript k and superscript N on the 
coefficients a{j;n). The time dependent Schrodinger equation for u{x,t) can be rewritten as a set of equations for 
these coefficients: 

A ^ 

[nuj + e- E°]a{j; n) ^ - E(i 1^1^) [«(^; n + 1) + a{l- n - 1)], (A8) 

1=1 

for j = 1,2, ... N. We emphasize that this is the exact equation for the time dependent problem in the finite basis 
N. It contains all orders of A and makes no reference to convergence of perturbation theory; there may be arbitrary 
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resonances or near resonances of the time dependent Hamiltonian between the initial N states. We will draw only on 
the fact that normalized solutions, satisfying ( |A8| ), exist. From Eq. (A7) we see that we need the coefficients a{j;n) 
only for the large values of frequency index n « Afjy. The high lying energies of the unperturbed static Hamiltonian 
are approximately Ej w i^^o (where uq = (1 / 2 m) ( tt / 2a)'^) . Then M^u) + e = -E'^v+i ~ i-^ + l)^<^o, and the factor 
in square brackets on the left hand side of Eq. ( |A8| ) for the case of interest, n — Mn, is greater than 2Nujq for any 
value of j (the smallest value occurs for the largest possible index, namely j — N). We use the symbol Vb to denote 
the maximum absolute value of the matrix elements of V between any two basis states. Then the absolute value of 
the coefficient a{j;n) is limited by Eq. ( [A^ ) to: 



Hj;Mn)\ < 



XVo 



N 



J2Mi;MN + l)\ + \aii;MN-l)\] < 



XVn 



(A9) 



where we have used only the limitation imposed by normalization, Eq. (A5), on sums over the absolute values of any 
subset of the coefficients a{j,n) corresponding to a single state , nam ely |ai| + \a2\ + ■ ■ ■ \am,\ < \An- But we can 
do much better, essentially by iterating this process. We start with Eq. (A8) for a smaller value of the photon index, 
n = Mn — p, use the argument just given to limit the right hand side for the next higher value, n — M„ — p + 1, 
and work back to the value of interest, n = M^. We choose the starting integer p as the integer part of Nujq/uj (this 
gives a substantial improvement only for Nloq/oj 3> 1 , so we choose N large enough for this to be the case). Then 
the coefficient in square brackets on the left side of (A8) is greater than iVwo- Thus, by exactly the same kind of 
argument that led to (A9) we have 



\a{3-MN-p)\ < 



(AlO) 



The same limitation holds for \a{j; Mn —p + 2q)\, with q = 1, 2, . . . ,p, where the coefficient on the left of (A8) is even 
larger. Then we use these maximum values to bound the right hand side of (AS) for the next iteration, for the values 
of n lying between those just limited: n = Mn — p + 2q — I with q = 1,2, ... ,p {not for q = 0): 



\a{j;MN-p + q)\ < 



XVo 
2NujQ 



N 



J2 MN-p + q+l)\ + \a{i; MN-p + q- 1) 



< 



XVo 



UJQ 



(All) 



where the first inequality comes again directly from (AS) with the minimum possible coefficient on the left hand side, 
and the final inequality results from substitution of (AlO), which holds for each value of the state index i, into the 
middle expression. This is repeated, using these bounds to limit the right hand side of (A8) for the n values lying 
between these: 



\a{j;MN -p + 2q)\ < 



1 



XVo 



1 3 



LUo 



(A12) 



now for q = 1,2, . . . ,p — 1. We then repeat this p — 2 times, with the power of (AVb/wo) increasing by one and the 
range of q decreasing by one at each iteration, to obtain 



Hj;Mn)\ < ^ 



XVo 



ujo 



P+i 



< 



1 



XVo 



ujo 



Nujo jo. 



Finally, we use this in Eq. (A7) to put limits on the size of the gap: 



A < AFoV2]V 



XVo 



0^0 



(A13) 



(A14) 



Therefore, the gap is limited by this to be exponentially decreasing with basis size iV, at least for small enough 
coupHng, A < wo/Vq. 

For specific examples of the spatial dependence of the time dependent potential V(x) cos ujt we can construct even 
tighter limits. There are, in particular, two limiting cases of interest: with the square well confining potential in the 
interval —a < x < awe take (i) V{x) = 2Vo sin(7rx/2a), or (ii) V{x) = VoaS{x). In both cases the only non-zero matrix 
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elements of V{x) between the eigenstates of the square well are of magnitude Vq. In the first instance (sinusoidal 
potential) these occur only for nearest neighbor states in the energy ladder, = VoS^ j^i, whereas for the delta 

function potential all pairs of even spatial parity states are connected by Vq regardless of how far apart in energy 
they are (and all odd parity states with vanishing wave function a,t x — 0, are, of course, totally unaffected by the 
potential). These are then limiting cases of short and long range effects of the spatial potential relative to the energy 

spectrum of the unperturbed static square well potential. 

For the sinusoidal potential t he ri ght hand side of Eq. (A8) contains only the four coefficients corresponding to 
I = j ±1, so that the inequality (A9) becomes 



|a(j;Mjv)| < 



2NuJo' 



(A15) 



The limitation to four coefficients occurs at each stage of the iterative process which led to (A13), which limit now 
becomes 



Hj;Mn)\ < ^ 



Nujo 



NujQ I u. 



(A16) 



Note that by choosing N large enough {Wq/Nujo < 1, as well as Nujq/lu > 1) we find ultimate exponential (indeed, 
powers of (1/7V!)) decrease of the gaps for arbitrarily large coupling strength A for this case. 



We also can obtain tighter limits for the delta function potential. In this case we can rewrite Eq. (A8) as 

m/2 



mo - 



[AN{n+l) + AN{n-l)], 



where we have defined 



N 



(A17) 



(A18) 



and throughout the analysis of this case the unperturbed eigenstate index j refers only to ev en pa rity states (we have 
noted above that the odd parity states arc not affected by this potential). Then we can sum ( A17 ) over the eigenstate 
index j to find a recursive relationship for the Aj^ (n) : 



Ajv(ri) = {XVo/2ujo)SN{n)[AN{n + 1) + Ajv(n - 1)], 
where we have defined one more sum. 



N r 



SN{n) = 



J 



N 

E 



1 



1 , Cn 

In ■ 



(A19) 



(A20) 



The final approximation on the right hand side is the Euler-Maclaurin integral estimate for the sum; corrections 
are of order 1/C„. Now, as before, we start by considering the recursion relation ( A19 ) for n — Mn — p, with p 
the integer part of Nljq /ld and limit the right hand side by the maximum imposed by the normalization condition: 

|^7v("-)| < V^, so that 



\An{Mn-p)\ < 



XVq \ii{AN + 1) 



/2N 



We use this in the right hand side of (A19) for the next higher value of n, namely n 
before to obtain 



(A21) 

Mn — p + I, and iterate as 



\An{Mn)\ < 



AFoln(4A^+ 1) 
2ujaV2N 



NujQ /a 



Finally, we put this back into the equation (A17) for the original coefficient a(j; Mjy) to find 



|a(j;A/jv)| < 



{2N + l)u;o 



XVa ln(4Ar + 1) 



- Nluq/lu 



2luoV2N 



(A22) 



(A23) 
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FIG. 1. Quasienergy spectrum as a function of A for A*' = 10 for the free particle in a box with harmonic driving and 
frequency to = 8.3 (see Eq. (5)). One finds many avoided crossings, with a typical one marked by the dashed box. There 
are some strictly real crossings, corresponding to states of opposite parity under the combined symmetry operation of spatial 
inversion plus time translation by half a period T/2. Other apparently real crossings are just so weakly avoided that they can't 
be resolved. 

FIG. 2. Position of branch points in the complex A plane for system and parameters of Fig. 1. Just one quadrant is shown, 
as the position of branch points is symmetric with respect to the axes Re(A) = and Im(A) = 0, since e(A*) — e*(A) and 
e(— A) — e(A). There are no branch points on the real axis, but they do appear on the imaginary axis, as shown. 
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